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Motivated by recent observational data, the equations of state with the inclusion of strangeness- 
bearing A-hyperons and the corresponding properties of neutron stars are studied, based on the 
covariant density functional (CDF) theory. To this end, we specifically employ the density de- 
pendent relativistic Hartree-Fock (DDRHF) theory and the relativistic mean field theory (RMF). 
The inclusion of A-hyperons in neutron stars shows substantial effects in softening the equation of 
state. Because of the extra suppression effect originated from the Fock channel, large reductions 



on both the star mass and radius are predicted by the DDRHF calculations. It is also found that 
t-H the mass-radius relations of neutron stars with A-hyperons determined by DDRHF with the PKA1 

parameter set are in fairly good agreement with the observational data where a relatively small neu- 
tron stars radius is required. Therefore, it is expected that the exotic degrees of freedom such as the 
strangeness-bearing structure may appear and play significant roles inside the neutron stars, which 
is supported further by the systematical investigations on the consistency between the maximum 
j^y^ neutron star mass and A-coupling strength. 

Q\ PACS numbers: 21.30.Fe, 21.60.Jz, 21.65.-f, 26.60.-c, 13.75. Ev 
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As the natural laboratories in the universe for nuclear 
and particle physics, neutron stars [T] have fascinated 
much effort concentrated on exploring the equation of 
state (EoS) of baryonic matter at low temperature and 
high density 13] . Specifically the mass of the observed 
neutron stars brings a strong constraint on the behavior 
of EoS at supranuclear density. The most precise mea- 
surements for the neutron star mass are determined to 
be less than 1.5M from the timing observations of radio 
binary pulsars [3] , which has remained for many years to 
constraint the EoS. However, the existence of more mas- 
sive compact stars is now unveiled by some evidence [5]. 
In a survey with the Arecibo telescope, an eccentric bi- 
nary millisecond pulsar PSR J1903+0327 was found with 
an unusually high mass value (1.74 ± 0.04)M Q [Hj. Re- 
cently a much larger pulsar mass of (1.97 ± 0.04)M Q was 
measured using Shapiro delay for the binary millisecond 
pulsar J1614-2230 [7]. All of these new data imply a 
stiff EoS of strongly interacting matter at high densities, 
which need further check on the new developed land- and 
space-based observatories. 

So far, there still exists considerable theoretical uncer- 
tainty on the EoS at supranuclear densities due to the 
poorly constrained many-body interaction, consequently 
deducing very different maximum mass and radius for a 
beta-stable neutron star. As indicated by prior study of 
neutron star based on the density dependent relativistic 
Hartree-Fock (DDRHF) theory jS], the maximum mass 
predicted by the covariant density functional (CDF) cal- 
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culations PHQ] lies between 2M and 2.8M . The Cor- 
responding EoSs deviate remarkably from each other at 
high density region. In the center of neutron stars, the 
density is generally considered as high as 5 to 10 times the 
nuclear equilibrium (saturation) density po ~ 0.16 fm~ 3 
of neutrons and protons found in laboratory nuclei. At 
such high density, exotica such as strangeness-bearing 
baryons, condensed mesons, or even deconfined quarks 
could come into existence [T] , which may play significant 
roles in determining the EoS. 

Besides the maximum mass limits, the mass-radius re- 
lation of neutron stars is also constrained by the recent 
observations, which leads to another strong restriction 
to the EoS. Recently, a relatively soft EoS and symme- 
try energy is predicted in the vicinity of the nuclear sat- 
uration density from careful analyses for three Type-I 
X-ray bursters with photospheric radius expansion and 
three transient low-mass X-ray binaries, which leads to 
a relatively small neutron star radius [TTJ [T2] . Such ob- 
servations are inconsistent with several commonly used 
equations of state that account only for nucleonic degrees 
of freedom. It is argued that they would be produced by 
including degrees of freedom beyond nucleons, e.g., hy- 
perons, mesons, and quarks, or possibly produced by a 
better description of nucleonic interactions [11]. As a 
possible solution without exotic degrees of freedom, the 
relativistic mean-field (RMF) models is recalibrated with 
a soft behavior of the symmetry energy around saturation 
density [HUH]. 

Compared to other CDF models, significant improve- 
ments have been obtained with DDRHF theory [TO] [15] 
in describing the relativistic symmetry conservation |15l - 
I17j . the consistency of the isospin dependence in nuclear 
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shell evolution [T51 Q]5] , the exotic structures [H] , and 
excited modes [22]. In the previous study [8], the signif- 
icant contributions to the symmetry energy have been 
found from the Fock terms in the isoscalar a and uj cou- 
plings, and the neutron star properties determined by 
DDRHF theory are shown in fairly good agreement with 
the observational data. 

In this paper, we continue the preivous work along the 
same line and study the roles of exotica in neutron stars 
based on the DDRHF theory. It is generally believed that 
hyperons appear around twice the normal nuclear matter 
density in neutron star matter 23 25j . The hrst hyperon 
to appear should be A as it is the lightest one with an 
attractive potential in nuclear matter. From the experi- 
mental binding energies of single-A hypernuclei, the po- 
tential depth of A in nuclear matter is estimated to be 
~ 30 MeV [26] ■ Motivated by the recent astrophysical 
observation, as a preliminary attempt it is thus interest- 
ing to introduce the strangeness degree of freedom asso- 
ciated with A-hyperon into the DDRHF theory and to 
investigate the corresponding neutron star properties. 

Within DDRHF, the baryons are described as point- 
like particles and interact with each other by exchang- 
ing mesons (the isoscalar a and uj as well as isovector 
p and 7r) and photons {A). The A-hyperon (A = uds), 
whose strangeness S = — 1, isospin 1 = and spin-parity 
J p = | + , participates only in the interactions propa- 
gated by the isoscalars, i.e., the a- and w-mesons. These 
interactions can be described by the following Lagrangian 
density, 

= -0a («7 M ^ - M A - g a . A a - g^.A^ui^) ip A , (1) 

where Ma denotes the mass of A-hyperon (^a)- From the 
Lagrangian density ([T]) the contributions to the energy 
functional, as well as the Dirac equation for A-hyperon, 
can be determined similarly as nucleon [5J [TUl [T5J H7J • 

For the beta-stable stellar matter containing nucleons, 
A-hyperons, electrons, and muons, the chemical equilib- 
rium conditions require that 

/Up =/i„ - /i e , (2) 

MA =P>n, (3) 

AV =Me (4) 

where the chemical potential p n , p,pi jj, A , fj,^, and p e are 
determined by the relativistic energy-momentum relation 
at Fermi momentum p — kp [5], 



Hi =E (k Fti ) + E*(k F ,i), H\=^k 2 FX + m\. (5) 

In above expressions, i denotes the baryons, i.e., neu- 
tron (n), proton (p), and A-hyperon, and A represents 
the leptons, electron (e - ) and muon (/z - ). Further com- 
bined with the baryon density conservation and charge 
neutrality, i.e., 

Pb =Pn + p P + Pa, P p =Ph + Pe, (6) 



the ratios of baryons and leptons can be obtained from 
the self-consistent calculations for the stellar matter with 
given baryon density p^. The equation of state (EoS) - 
the relation between the pressure and energy density for 
stellar matter, is then determined by DDRHF as well as 
other CDF calculations, in which the baryons interact by 
exchanging mesons whereas the leptons are described as 
free fermions. 

With the EoS of the stellar matter, the structure of a 
static, spherically symmetric, and relativistic star can be 
determined by solving the Tolman-Oppcnhcimcr-Volkov 
(TOV) equations [Ml [^.Similarly as in Ref. [5], the 
EoS in the low density region (p& < 0.08fm~ 3 ) will be 
provided by the BPS [3D] and BBP 31] models. For a 
given central density p(0) or central pressure -P(O), the 
input EoS leads to the unique solution of the TOV equa- 
tions. 

In the following the theoretical calculations are sys- 
tematically performed on the platform of CDF theory 
[lOj [111 EZl E2l [33], specifically the DDRHF theory with 
effective interactions PKAl [T5J and PK03 [T5] , and the 
RMF theory with NL-SH (SJ, PKl [35], TW99 [3J, and 
PKDD [35] . In Table [i] are summarized the details of 
the selected effective Lagrangians. For the RMF calcu- 
lations, only the Hartree contributions are involved in 
the baryon-baryon interactions, whereas both Hartree 
and Fock terms are taken into account by DDRHF. In 
NL-SH and PKl, the non-linear self-couplings of a- and 
w-mesons [35J [371 13H] are introduced to evaluate the in- 
medium effects of nuclear interaction. In TW99, PKDD, 
PKAl, and PK03, such effects are introduced by the 
density dependence in meson-nucleon couplings |39|, I40j 
and the same density dependent behaviors are utilized 
in corresponding meson-A coupling channels, i.e., the 
meson-baryon coupling constants are treated as functions 
of baryon density (pb = p n + Pp + Pa)- As proved by prior 
studies [5], such prescriptions on the in-medium effects 
are still reasonable for exploring the EoS at supranuclear 
density as well as in describing the neutron stars. 



TABLE I: Details of the selected CDF effective Lagrangians. 





Fock 


non-linear 


non-linear 


density dependent 




term 


a term 


w term 


couplings 


PKAl 


yes 


no 


no 


yes 


PK03 


yes 


no 


no 


yes 


NLSH 


no 


yes 


no 


no 


PKl 


no 


yes 


yes 


no 


TW99 


no 


no 


no 


yes 


PKDD 


no 


no 


no 


yes 



In the CDF calculations if not specified the propor- 
tions between meson-A and meson-nucleon couplings are 
fixed as g a -Ajga — 0.600 and g^-A/gu = 0.653 [41], 
and the masses of A-hyperon, electron, and muon are 
chosen as M A = 1115.0 MeV, m e = 0.511 MeV, and 
to m = 105.7 MeV. In the Lagrangian density ([!]) we ne- 
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FIG. 1: (Color online) The pressure as a function of the bary- 
onic density pb (fm -3 ) with different CDF effective interac- 
tions for the stellar matter containing nucleons, A hyperons, 
electrons and muons (red curves), as compared to the one 
without A hyperons (black curves). The corresponding contri- 
bution from Hartree and Fock channel with PKA1 are shown 
in the lower panel in comparison with TW99. See the text for 
details. 

gleet the strangeness concerned baryon-baryon interac- 
tions, which still remain to be an open question. In fact, 
several recent observations of double-A hypernuclei indi- 
cate that the effective AA interaction should be consid- 
erably weaker than earlier evaluations |42) . 

With the selected CDF effective Lagrangians, we first 
study the equation of state for the beta-stable stellar 
matter containing nucleons, A-hyperons, electrons and 
muons (or NAep for short), in comparison with those 
without A-hyperons (or Nep for short). Fig. [ija) shows 
the pressures of neutron star matter as functions of the 
baryon density p b with the selected CDF effective interac- 
tions. For the Nep matter, there exist substantial devia- 
tions (black shadowed ellipse) between the EoSs, among 
which NL-SH presents the hardest EoS and TW99 gives 
the softest one. When the A-hyperon degree of freedom 
is introduced, the EoSs become much softer above the 
critical density (~ 0.3 fm -3 ), where A-hyperons start to 
appear in stellar matter. In addition, the deviations (red 
shadowed ellipse) on the EoSs are considerably reduced 
as well. The distinctly stiff behavior of EoS with NL- 



SH effective interaction is mainly due to the absence of 
the non-linear ui self-couplings [37]. It is also found that 
PKA1 presents rather harder EoS than TW99 in the Nep 
matter, while with A-hyperons PKA1 and TW99 tend to 
provide almost identical EoS until p b ~ 0.8 fm -3 , where 
the deviation emerges again. 

In fact, such features on the EoS can be understood 
qualitatively from the chemical equilibrium between nu- 
cleons and hyperons [see Eq.([3|]. For harder EoS, the 
neutron chemical potential goes up faster with the den- 
sity increasing, such that more neutrons will be trans- 
ferred into A-hyperons. Because the w-meson couplings 
with A-hyperons, which play dominant role at high den- 
sity, are generally weaker than those with nucleons as re- 
pulsions, the appearance of A-hyperon will soften the EoS 
evidently. Thus, more stiffer a EoS of the Nep matter 
is, more effort to soften the EoS is made by A-hyperon. 
When the density reaches at rather large values, the ra- 
tios of A-hyperons and nucleons become stable so that the 
deviations among the effective Lagrangians appear again. 
Even though, such deviations have been diminished ev- 
idently due to the fairly large proportion of A-hyperon 
at high density. Compared to the constraint [shadowed 
area in Fig. [IJa)] suggested in Ref. [11] for cold matter, 
the EoSs of NAep matter show much better consistency 
with the constraint than those for Nep matter. 

To illustrate the influence of the Fock term on the EoS 
for the Nep and NAep matter, the pressures of neutron 
star matter with DDRHF effective interaction PKA1 and 
DDRMF one TW99 are replotted in Fig. [TJb) as func- 
tions of the baryon density. The contributions from the 
Hartree and Fock channels are shown as well. Within 
DDRHF, one can find that the Hartree terms provide 
the dominant contributions to the pressures in compar- 
ison with the Fock terms in both Nep and NAep mat- 
ter. It is also seen that the Hartree contributions in 
PKAl are nearly identical with TW99 and the differ- 
ence in the EoSs comes significantly from the Fock term 
in PKAl. Compared to the Nep matter, the occurrence 
of A-hyperons remarkably suppresses the Fock contribu- 
tions in the NAep matter. This is mainly due to the 
effects of uj couplings in the Fock channel. In Nep mat- 
ter the Fock terms in the N-ui couplings contribute as 
fairly strong repulsion at high density. Whereas in NAep 
matter the A-u> couplings in Fock channel represent as at- 
traction, which substantially softens the EoS. As a result 
nearly identical behavior of EoSs is provided by PKAl 
and TW99 in the density region about 2 ~ 5po, showing 
better consistency with the constraint [TT] than other 
CDF effective Lagrangians. 

The symmetry energy is an important quantity for 
illustrating the property of asymmetric nuclear mat- 
ter. In general, the energy per particle of asymmet- 
ric nuclear matter, E(p b ,f3), can be expanded in the 
Taylor series with respect to the asymmetry parameter 

P= (Pn~ Pp)/{Pn + Pp), 

E{p b , p) = E { Pb ) + f3 2 E sym (p b ) + ■■■ . (7) 
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FIG. 2: (Color online) The symmetry energy E aym (MeV) as 
a function of the baryonic density pb (fm~ 3 ) with different 
CDF effective interactions for the stellar matter containing 
nucleons, A hyperons, electrons and muons (red curves), as 
compared to the one without A hyperons (black curves). The 
corresponding contribution from Hartree and Fock channel 
with PKA1 are shown in the lower panel in comparison with 
TW99. See the text for details. 



The function E$(p b ) is the binding energy per particle in 
symmetric nuclear matter. The empirical parabolic law 
in Eq. Q is confirmed to be reasonable throughout the 
range of the asymmetry parameter values especially at 
low density. As a reasonable approximation, one can ex- 
tract the symmetry energy E sym (pt,) for the beta-stable 
stellar matter by 



E sy m{pb) — 



E(p b ,P)-E (p b ) 
P 2 



(8) 



Fig. [2ja) shows the symmetry energy of neutron star 
matter as functions of the baryon density pt> with dif- 
ferent CDF effective interactions. For the Nep matter, 
sizable enhancements of E sym are seen in the high den- 
sity region with respect to E sym at the saturation denstiy. 
When the A-hyperon is introduced, the symmetry ener- 
gies become much softer, about 50% reduced at high den- 
sities. Similar as the results in Fig. [T] the deviations of 
E sym among different EoSs are also reduced in the NAep 
matter. For the Nep matter, the Fock terms of a- and lo- 
couplings exhibit significant contributions to the symme- 



try energy, which leads to a stronger density dependence 
in DDRHF than in DDRMF at high density [5]. While 
in the NAep matter, both Hartree and Fock contribu- 
tions (e.g., in PKA1) are much reduced by A-hyperons. 
In Fig. [2jb) one may find that the Hartree contribu- 
tions in PKAl are essentially identical with TW99 and 
correspondingly the Fock terms provide about 2-3 times 
reduction as the Hartree ones with the occurrence of A- 
hyperons, leading to notably soft symmetry energy at 
high densities. 

Taking the EoSs in Fig. [lja) as the input, the sys- 
tematic properties of neutron stars such as the mass- 
radius relation can be obtained by solving the TOV equa- 
tions. In Fig. [3]ja) are shown the mass-radius relations 
for the neutron stars which consist of the NAep mat- 
ter (solid lines), compared to the Nep matter (dashed 
lines). Consistently with the EoSs in Fig. [I] the in- 
clusion of A-hyperon brings substantial reduction on the 
neutron star mass. When A-hyperons appear in neutron 
star, the mass-radius relations deviate from those with- 
out A-hyperons and bend down to smaller masses and 
radii. Close to the maximum points (in symbols) the 
mass-radius relations tend stable, which is also different 
from the calculations without A-hypcron. This can be 
interpreted by the behaviors of symmetry energy at high 
density. In Nep matter the symmetry energies increase 
with the density almost linearly. While with the inclusion 
of A-hyperons the symmetry energies bend down to be 
stable or even slightly decreasing when the density goes 
high, in a consistency with the mass-radius relations [see 
Fig. §a)]. 

Moreover, the CDF calculations for the Nep matter 
predict quite different mass for neutron stars. How- 
ever, when A-hyperon is included, the deviations of the 
maximum masses among different EoSs are consider- 
ably reduced, namely, the CDF effective Lagrangians ex- 
cept NL-SH predict vicinal maximum mass. Comparing 
PK03 with PKDD and PK1, one can find that nearly 
identical maximum masses are provided by PKDD and 
PK1 in both cases whereas PK03 gives larger value in 
the case of the Nep matter. As mentioned above this can 
be understood well from the effects of the Fock terms in 
softening the EoS and thus reducing the maximum mass 
of neutron stars. 

Extracted from Fig. [3]ja), Table [II] shows the mass 
(Mq), radius (km), and central density p c (fm -3 ) for neu- 
tron stars with the maximum mass, as well as the thresh- 
old densities p b (fm -3 ) and p b (fm -3 ) for A and muon 
emergence. For comparison, the results given by the 
calculations without A-hyperon are shown in the lower 
panel. From Table [H] one can clearly see the mass reduc- 
tion induced by the occurrence of A-hyperons, especially 
in the CDF calculations with the Fock terms, i.e., the 
DDRHF calculations (PKAl and PK03), which present 
the mass reduction about 0.7M Q . In the RMF calcu- 
lations the mass reductions range from 0.4 to 0.5 Mq, 
except NL-SH which gives the reduction about 0.6M©. 

The inclusion of A-hyperon also brings distinct ef- 
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FIG. 3: (Color online) The mass-radius relations for the 
neutron stars with A (solid lines) and without A (dashed 
lines) hyperons. The symbols denote the neutron stars with 
maximum masses. Region excluded by causality is indi- 
cated with gray color. For comparison the observational con- 
straints from three neutron stars in the binaries 4U 1608—248 
(green/red), EXO 1745-248 (yellow/blue), and 4U 1820-30 
(cyan/ magenta) in Ref. [1 11 are shown as the 1- and 2-cr 
confidence contours, and the later analyses in Ref. |12| are 
denoted by the shadowed areas with two different models of 
photospheric radius. In the lower panel are shown the mass 
& radius for the neutron stars with maximum masses deter- 
mined with different A hyperon coupling strength (filled sym- 
bols), as compared to the ones without A hyperons (open 
symbols). See the text for details. 



fects on central density p c . Among the selected effec- 
tive Lagrangians PKA1 and TW99 predict the largest 
values of p c , about 8 times of the saturation density (~ 
0.16 fm ). Compared to those excluding A-hyperon, 
the central densities increase roughly 25%~38% (about 
0.25 ~ 0.50fm -3 ) in the DDRHF calculations (PK03 and 
PKAl), and about 10%~15% in the density dependent 
RMF calculations (PKDD and TW99). While the central 
density predicted by PK1 decreases with the A-hyperon 



TABLE II: Maximum mass M max (Mq), corresponding radii 
iJmai (km) and central density p c (fm -3 ) for neutron stars, as 
well as the threshold densities p£ (fm -3 ) and (fm -3 ) of A 
and muon. For comparison, the quantities for neutron stars 
without A hyperons are shown in the lower panel. The results 
are calculated by DDRHF with PKAl and PK03, and RMF 
with PKDD, TW99, PK1, and NL-SH. 





PKAl PK03 


PKDD TW99 


PK1 NL-SH 


^max 
^max 
Pc 


1.713 1.837 
10.425 11.495 
1.314 1.048 


1.849 1.647 
11.583 10.333 
1.024 1.292 


1.832 2.213 
13.048 13.633 
0.786 0.700 


pt 
Pt 


0.272 0.284 
0.118 0.122 


0.322 0.368 
0.108 0.116 


0.306 0.282 
0.110 0.114 


■^max 
Pc 


2.423 2.500 
12.354 12.487 
0.810 0.780 


2.329 2.078 
11.798 10.632 
0.888 1.104 


2.315 2.802 
12.705 13.534 
0.796 0.650 



inclusion and NL-SH simply gives tiny changes in p c . 

For the radii of neutron stars with the maximum mass, 
distinct reductions are also found in the DDRHF calcula- 
tions including A-hyperon, e.g., PKAl predicts the star 
radius about 10.4 km, about 2 km smaller than the one 
without A-hyperon, and PK03 provides the reduction 
about 1 km. In the RMF calculations, the neutron star 
radii are changed only slightly (< 0.3km) with the in- 
clusion of A-hyperon. Such discrepant behavior between 
DDRHF and RMF originates from the correlations be- 
tween the radius of neutron stars and the symmetry en- 
ergy. In Table |III| are shown the bulk quantities of sym- 
metric nuclear matter at saturation density for the se- 
lected effective Lagrangians. As seen from the slope L 
and curvature K sym on the symmetry energy [15], the 
DDRMF calculations (PKDD and TW99) present rather 
soft symmetry energy and the slopes tend to decrease as 
density goes high due to the negative curvature K sym . 
In an accordance, small neutron star radii are obtained 
in the case of Nep matter, especially for TW99. Simi- 
lar consistence can also be seen from the calculations of 
DDRHF and NLRMF. In Fig. [2] it is clearly shown that 
the occurrence of A-hyperons brings substantial effects 
in reducing the slope of the symmetry energy, mainly 
from the Fock channel in DDRHF. Due to such extra 
suppression, consistently the reductions of the star ra- 
dius are more dramatic than the case of RMF, as proved 
by the results in Table [H] On the other hand it is also 
well demonstrated that the properties of neutron star 
are strongly correlated with the values of the symme- 
try energy. As seeing from the results in Fig. [2ja), one 
should notice that the symmetry energies are dramati- 
cally changed with the occurrence of A-hyperon at high 
density so that the properties of symmetry energy at nor- 
mal density are not enough to describe properly the size 
and mass of neutron stars. 

In Fig. |3j the 1- and 2-cr confidence contours for the 
masses and radii of three neutron stars in the binaries 
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TABLE III: Bulk properties of symmetric nuclear matter 
at saturation point, i.e., the saturation density po (fhi~ ), 
binding energy per particle E/A (MeV), incompressibility K 
(MeV), and symmetry energy J (MeV) with its slope L (MeV) 
and curvature K SYm . (MeV). The results are provided by the 
CDF effective Lagrangians PKA1, PK03, PKDD, TW99, NL- 
SH, and PK1. 





Po E/A K 


J L ^sym. 


PKA1 
PK03 


0.160 -15.83 229.96 
0.153 -16.04 262.47 


36.02 103.50 213.23 
32.99 83.00 116.56 


PKDD 

TW99 


0.150 -16.27 262.18 
0.153 -16.25 240.27 


36.79 90.21 -80.74 
32.77 55.31 -124.68 


NL-SH 
PK1 


0.146 -16.35 355.43 
0.148 -16.27 282.69 


36.12 113.66 79.72 
37.64 115.88 55.33 



4U 1608-248 (green/red), EXO 1745-248 (yellow/blue), 
and 4U 1820—30 (cyan/magenta) are shown [TTJ, and a 
re-analyzed version recently done by another group is 
denoted by the shadowed areas with two different mod- 
els of photospheric radius as well [12]. Different from 
other constraints like shown in Ref. these new ob- 
servations put a serious challenge to our understanding 
of neutron stars, where a relatively small values about 
8.7-12.5 km is required for a 1.4M Q star, even smaller 
than a recent conclusion of 9.7-13.9 km from microscopic 
calculations based on chiral effective field theory interac- 
tions [44] . Therefore a very soft symmetry energy and 
EoS near and above several times of the saturation den- 
sity are needed. In comparison with our calculations, it 
is found that the mass-radius relations for the neutron 
stars with A-hyperons given by PKA1 and TW99 are in 
perfect accordance with the observations, especially for 
the cases around the maximum mass. Because of a little 
harder EoSs at the density region about 2 ~ 5po, PK03 
and PKDD with A-hyperons just marginally cover the 
constraints, while the nonlinear RMF versions PK1 and 
NL-SH could not fulfill the constraints at all. In the cases 
without A-hyperons, all the curves are far away from the 
constraints. Hence, it is strongly suggested the exotic 
degrees of freedom such as the strangeness-bearing struc- 
ture may appear inside the neutron stars. It is expected 
and also found in the CDF calculations that A-hyperon 
is the dominant constitution in the core of neutron star 
whereas neutron is strongly compressed to be less than 
20%, from which the role of the hyperon degree of free- 
dom is well demonstrated in neutron stars. 

In Fig. [I] and Fig. [3]Ja) it is already shown that the 
A-hyperon plays an important role in softening the EoS 
and reducing the neutron star masses and radii and leads 
to fairly good agreements with the constraints [HI H2] • 
One may notice that in the case of NAe/i matter none 
of the curves except NL-SH go through 2.0M Q , which 
is constrained by another observation [7j. On the other 
hand if comparing the results provided by the selected 
effective Lagrangians it is also found that for the CDF 
calculations with Fock terms, i.e., the DDRHF ones, the 



inclusion of A-hyperon brings more distinct effects in re- 
ducing the neutron star masses. 

In all the above calculations the coupling strengths of 
A-hyperon are fixed to <? ct -a/<?ct = 0.600 and g u -hl9u> = 
0.653. For the A-hyperon, it only participates in the in- 
teraction mediated by the exchange of the isoscalar <j- 
and w-mesons, which respectively represent as strong at- 
traction (repulsion) and repulsion (attraction) in Hartree 
(Fock) channels. It is well known that at high density the 
contributions from w-meson play the dominant role in de- 
termining the EoS, as well as the mass-radius relation for 
neutron stars. In DDRHF there exist substantial cance- 
lations between the Hartree and Fock terms of the w-A 
coupling and such cancelations are somewhat equivalent 
to weakening the coupling strength. More distinct effects 
are therefore found in the DDRHF calculations with the 
inclusion of A-hyperon in softening the EoS, as well as 
reducing the neutron star mass. 

Here we only consider the strangeness system degree 
of freedom associated with A-hyperon. Approximately 
and qualitatively, effectively other strangeness related 
contributions such as E- and S-hyperons can be taken 
into account by modifying the A-coupling strength. In 
Fig. |3jb) are shown the neutron stars with the max- 
imum mass extracted from the CDF calculation with 
different A-coupling strengths, where the Greek symbol 
A denotes the coupling strengths g a -A/g<7 — 0.600 and 
Su-A/ffw = 0.653. That is, for 1.5A for instance, the cou- 
pling strengths are taken to be g a -\/ga = 1-5 x 0.600 
and g^-h/guj — 1-5 x 0.653. It is found that the neutron 
star masses are substantially reduced as the A-couplings 
change from 1.5A to 0.5A, roughly corresponding to the 
uncertainty in A-coupling. Such behavior can be in- 
terpreted by the consistency between the EoS and A- 
coupling strength. With the weakening of A-coupling, 
which is equivalent to reducing the A repulsion as well as 
the Fermi energy, more and more neutrons will be trans- 
ferred into A-hyperons and the EoS will become softer 
and softer. From Fig. [3jb) one may find that the 2M 
constraint is fulfilled reasonably despite the uncertainty 
of A-coupling. It is also found that with reducing A- 
coupling the radii of neutron stars decrease to the min- 
imum first and then keep increasing. The CDF calcu- 
lations also show that the central densities increase to 
the maximum first and then keep decreasing with the re- 
duction of A-coupling, being consistent with the radius 
evolutions. As shown in the plot such turning points in 
the mass-radius relations are rather close to the original 
ones (denoted by A), especially for PKAl, which may 
imply that the coupling strengths ga-A/ga — 0.600 and 
5u-a/.9w = 0.653 are reasonable for the A-coupling in stel- 
lar matter as well as in finite nuclei |41j . 

In summary, we studied the general properties of 
neutron stars with the inclusion of strangeness-bearing 
A-hyperon, based on the covariant density functional 
(CDF) theory, specifically the density dependent rela- 
tivistic Hartree-Fock theory and the relativistic mean 
field theory with both nonlinear self-coupling of mesons 
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and density dependent meson-nucleon couplings. The in- 
clusion of A-hyperon in neutron star systems shows sub- 
stantial effects in softening the equation of state for the 
stellar matter, as well as in reducing the star mass and 
radius, especially when the contribution of the Fock chan- 
nel is included. It is shown that the properties of symme- 
try energy at normal density are not enough to predict 
the mass and radius of neutron star. The systematical in- 
vestigations on the consistence of the maximum neutron 
star mass and A-coupling strength also indicate that ex- 
otic degrees of freedom is one of the important factors to 
provide appropriate prediction on the star mass as well 
as the radius in consistency with recent observations. 
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